Coherent tunneling and quantum coherence oscillations 
at the atomic level 1 
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The evolution of the quantum wave packet describing an atom trapped in the 
surface-tip junction of the scanning tunneling microscope is investigated by 
using the time-dependent Schrodinger equation, and a quasi-classical Hamil- 
tonian approach. The estimates concern a Xe atom in a biased double-well 
junction potential. The exact treatment shows that quantum coherence os- 
cillations of the metastable ground state may occur at particular resonant 
values of the bias voltage. The effect of decoherence by partial localization 
is studied within the quasi-classical frame. 
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1 Introduction 



The phenomenology of the wave function collapse and decoherence is a sub- 
ject of wide interest, ranging from the conceptual framework of quantum 
mechanics [1, 2] to the physics of the quantum logic circuits. 

Dynamical loss of coherence produced by possible non-linear terms in the 
time-dependent Schrodinger equation (TDSE) could be observed by measure- 
ments on driven hyperfme transitions in 9 Be + ions [3]. Such non- linear terms 
can appear by the coupling to the environment [4], but they could also be in- 
trinsically built into the quantum mechanics, as it was proposed by Ghirardi, 
Rimini and Weber to explain the paradoxes of the measurement theory [5]. 
The mean-field of molecules or atomic nuclei provided by the Hartree-Fock 
approximation breaks spontaneously the translation symmetry, and is local- 
ized in the inertial frame defined by the center of mass (CM). The CM motion 
can be quantized, removing the spurious effects of the symmetry breaking, 
but when the size of the objects increases towards the macroscopic level, (the 
object becomes "environment" for its constituents), physical localization may 
occur [6, 7]. 

The interplay between quantum coherence and wave function collapse ap- 
pears particularly striking in the computation models of the quantum com- 
puters [8]. Therefore, a most interesting issue is to study the coherence 
properties of the physical systems at atomic level. Recent experiments have 
proved the existence of a coherent component for the electron tunneling in 
quantum dots [9], while quantum localization was predicted for atoms mov- 
ing in the periodic potential created by a phase-modulated standing light 
field [10]. 

The purpuse of this work is to investigate the feasibility of experiments 
on localization and decoherence for the CM wave function of the individual 
atoms, using the scanning tunneling microscope (STM). Since the first exper- 
iments on reversible atom transfer [11], STM may be considered as the ideal 
instrument for manipulating atoms or molecules. The bistable operation 
mode of STM can be understood assuming that the diffusion barrier on the 
surface is high enough to prevent the escape of the particle from the junction 
region, and that the motion takes place along the outer normal to the surface 
plane in an asymmetric, one-dimensional, double- well potential (DWP) [12]. 
In a DWP, an observable effect which is very sensitive to decoherence is the 
tunneling of a wave-packet tpo with the energy below barrier [13]. If ipo is a 
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Gaussian localized initially in the metastable well, and the quantum coher- 
ence is preserved in time, then complete tunneling to the stable well appears 
only in very special resonance conditions, by quantum coherence oscillations 
(QCO) [14]. 

The coupling to the environment affects the superposition principle and 
the quantum coherence property, and in the case of a dissipative force linear 
in the CM velocity, the QCO between the two wells may either be damped 
[15], or completely suppressed [16]. However, decoherence could be produced 
also by dissipative terms which only suppress the spreading of the wave packet 
(e.g. the "squeezing" dissipation [4]). In the asymmetric DWP, such terms 
may change the tunneling mechanism from QCO between wave functions 
with different shapes, to coherent tunneling (i.e. after barrier crossing the 
wave function retains its initial shape [14]). 

In Sect. 2 it is shown that for a Xe atom placed in the STM surface-tip 
junction, reversible tunneling by QCO resonances could be observed. The 
irreversible switching mode is discussed in Sect. 3, within a quasi-classical 
model. Conclusions are summarized in Sect. 4. 

2 The atomic quantum coherence oscillations 

An isomeric state tp of a quantum particle in an asymmetric DWP has 
a QCO resonance if it is a linear superposition of two quasi-degenerate 
eigenstates ipd, ipu of the Hamiltonian, with energies ~ E u , such that 
A = \E U — Ed\ is much smaller than the average level spacing. In this case 
the localization probability in the stable well, p(t) , oscillates according to the 
law p(t) = [1 — cos(irt/T max )]/2, with T max = Tin/ A. In most physical situ- 
ations the resonance conditions are not fulfilled, but particularly interesting 
are the systems where the potential depends on a parameter which may be 
changed continuously until there is resonance 2 . For STM this parameter can 
be the surface-tip bias voltage U. 

The present estimates are based on the potential function given in ref. 
[12] for Xe atoms in a fixed surface-tip geometry. This potential consists in 
a DWP term produced by the binding interaction energy, V p (x), (Fig. 1, 
dashed line) and the dipole term, Va(x) = — U/i {[0.3 + 0.7(w + x) 4 / ' L 4 }' 1 — 

2 Resonances in a three- well potential of variable width are presented in the application 
2349433, Long distance quantum transfer, available at http://cipo.gc.ca 
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[0.3 + 0.7(w - x) 4 /L 4 ]~ 1 }/2w, where /i = 0.3 D, (lD=3.335x 10~ 30 C-m), 
w = 2.2 A, and L = 1.56 A. The X-axis is chosen normal to the surface, with 
the origin at the barrier top of the term V p . This term is taken as a fourth or- 
der polynomial, with the metastable minimum Vq — Vp(%o) — —12.7 meV at 
Xq = —0.7 A, near surface, the barrier top V& = V p (xb) = 0.45 meV at xj, = 0, 
and the stable minimum V g = V p (x g ) = —23.2 meV at x g = 0.89 A, close 
to the tip. At small bias the total potential V(x) = V p (x) + Va(x) remains 
a DWP, with the coordinates xo,Xb,x g of the extrema close to the values 
given above. For negative bias of the surface the barrier height — Vb — Vq 
decreases, and at U ~ -1.2 V the isomeric minimum dissapears. 

The CM wave function for the metastable ground state of a Xe atom 
localized initially near surface is well approximated by the Gaussian wave 
packet ipo{%) = (c /7r) 1//4 e~ Co ^ _:Eo ) / 2 , where c = Mu /h, M is the Xe mass, 
and u = ^V"{xq)/M is the classical oscillation frequency at x in the har- 
monic approximation. 

If there are no decoherence factors, the wave function at the moment t is 

1,(x,t) = e-rt^M*) , *i = -^^ + V p (x) + V d (x) , (1) 
and can be obtained by integrating the TDSE 

ih^ = m (2) 

with ip(x, 0) = ipo(x) as initial condition. 

The TDSE can be solved numerically in a spatial grid {xk}, k=l,N, using 
the leap-frog method [17]. However, the computation becomes much faster if 
the discrete form of Eq. (2) is integrated as a Hamilton system of equations. 
Thus, if Uk(t) = Re(ip(xk,t)) and Vk(t) = Im(ifj(xk,t)) denote the real, 
respectively the imaginary part of the wave function if)(x, t) at the grid point 
Xk, then Eq. (2) becomes 

9n ^. &h 

2hu k = — 2hv k = -- — , (3) 

ov k du k 

with 

N 

H = ^2u k (Hu) k + v k (Hv) k , (4) 
fe=i 
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W)* = I 90 3 20 + 3 2 T8 yfc] + 

+[V p (x k ) + V d (x k )\y k . 

In the present application the Hamiltonian system of Eq. (3) was defined 
considering N = 321 spatial grid points equally spaced by I = 0.01 A within 
the interval [x m i n , x max \— [—1.2 A, 2 A]. For a fast integration can be used 
the D02BAF routine of the NAG library [18], with the time step dt = 
6.58 x 10~ 2 ps. The localization probability in the stable well of V(x) is 
P(t) = fx™* ^*{ x i t)tp(x,t)dx, and within the time interval [0,20 ps] attains 
a maximum which is represented in Fig. 2(A) as a function of the bias 
voltage U. The peaks indicate the QCO resonances. The first appears at 
U = —1.141 V for the potential V represented in Fig. 1 by solid line. This 
resonance corresponds to p(t) shown in Fig. 2(B), with the maximum at the 
moment T max = 14.37 ps. In Fig. 1 can be seen the wave function at t — 0, 
(ip , dotted line), and at t — T max , (ipM, solid line). 



3 The coherent tunneling 

During QCO the wave-packet changes its shape, and therefore the tunnel- 
ing is not coherent. The issue of coherent tunneling may receive an answer 
from the study of the evolution pattern for wave functions constrained to 
be Gaussian all the time, as in the quantum molecular dynamics [19]. The 
constrained dynamics will be obtained by a time-dependent variational cal- 
culation within the trial manifold of the Gaussian wave packets with variable 
centroid and width. This manifold contains the isomeric ground state, and 
is represented by a combination between coherent and squeezed states, 

i> c {z, s) = e (^^b) e ( sb W-s^ b )/2^ o ? z = a + ifr s = pe-** , (5) 



with (a, /3, p,(f>) real parameters, and b = y Mujq/2Ti(x — x + hd x /Muo) 
the Dirac-Fock annihilation operator for ip . In terms of the variables x c = 
xo + a^2h/Muo, p c = (3\/2hMuo, v = (sinhp) 2 , the equations of motion 
obtained from the variational equation 5 J dt(ip c \ihd t — H\ip c ) = are 

d(ij c \H\ij c ) . _ d{^ c \H\^ c ) 



Xr , , ; 

op c dx c 
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• d(^ c \H\^ c ) d(^ c \H\jj c ) 
H =—dv hv = 90—- (7) 

Eqs. (6) describe the motion of the CM coordinate and momentum x c , p c of 
the Xe atom, while Eqs. (7) provide the evolution of the localization width 
around x c , a 2 = (K) = (x 2 ) - (x) 2 = [1/2 + v + Jv(v + l)cos20]/c o . To 



express the average (ip c \H\ip c ) in analytical form, the potential V(x) was re- 
placed by a fourth order interpolation polynomial, V po i(x), with coefficients 
determined by fit. 

At the first resonance (U = —1.141 V), the best fit of the potential V(x) 
(Fig. 1, solid line) is given by a quartic polynomial V po i(x) (Fig. 3(A), 
solid line), with the extremum values of 0.8 meV, 1.85 meV and —48.88 
meV at the points x = —0.47 A, xt = —0.19 A, and x g = 1.01 A, respec- 
tively. The classical oscillation frequency at xq which defines ipo in Eq. (5) is 
c^o — 1-95 ps~ x . Using these results, H at resonance is well approximated by 
Hp i = -h 2 dl/2M + V po i(x), and V av (x) = (H po i)(x c = x,p c = 0,v,(f>) has the 
role of effective potential for the CM dynamics of the Gaussian wave packet 
ip c . This function is represented for = and v — 0, 0.1, 0.5, 1, 1.5, 2 in Fig. 
3(A) by dashed lines. If v = 0, then V av (x) is also a DWP, with the extremum 
values of 1.57 meV, 1.97 meV and -46.8 meV at -0.425 A, -0.215 A, and 
1 A, respectively. Thus, for a Gaussian wave packet the effective potential 
is not the same as for the classical particle, and the top of the barrier has a 
small shift to higher energy. 

When x c = xq and p c = v = <p — 0, the trial function ip c reduces to 
the isomeric ground state i/jq of H po i- The quasi-classical dynamics of a Xe 
atom which is initially in this state was obtained by integrating the system 
of Eqs. (6), (7) using the routine D02BAF of the NAG library, and the tra- 
jectory x c (t) is represented in Fig. 3(B) by dashed line. For comparison, the 
exact wave function ip{t) = exp(—itH po i / h)tp was calculated by solving the 
TDSE, and the corresponding expectation value XQcoif) — (V ; WM' ? /'W) is 
represented in Fig. 3(B) by solid line. Initially x c (t) and xqcoif) ar e close, 
but in time XQcoif) extends across the barrier, while x c (t) oscillates in the 
isomeric well with small amplitude 3 . 



3 Numerical estimates indicate that if the "squeezing" dissipation term provided by the 
coupling operator K = (x — {xp\x\ip)) 2 is included, and ip(t) is obtained by integrating the 
non-linear TDSE ihdtip — (H po i + nKdt(ip\K\ip))i[>> then for k ~ 50h/ A 4 any significant 
changes in the shape off/' are suppressed, and (ip(t)\x\ip(t)) remains close to x c {t). 
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The oscillation amplitude of x c increases with the energy E c = (ip c \Hpoi\ipc) , 
but the trajectory remains confined in the isomeric well until E c reaches a nar- 
row interval very close to the barrier top of V av , when a "switching" mode ap- 
pears. For an orbit with the initial conditions x c ~ —0.25 A, p c = v = <fi = 0, 
the oscillations in the isomeric well are interrupted by a sudden jump to the 
stable well, where the particle keeps oscillating without return. This inter- 
esting behaviour resembles the experimental situation, because the observed 
atom switching occurs with a rate proportional to a power of the current, 
and could be explained by the electron heating effects [11, 12]. Here, this 
mode is due to a small valley connecting the two wells which appears near the 
barrier top of V av in the squeezing dimension. Thus, the escape to the stable 
well becomes possible when the four dimensional orbit C t = (x c ,p c , v, </>)(£) is 
directed along this valley. The projection of the switching mode orbit C on 
the CM phase space plane is pictured in Fig. 4(A), while the corresponding 
trajectory plot, x c (t), is shown in Fig. 4(B). 

4 Conclusions 

The TDSE calculations indicate that if quantum coherence is preserved, then 
elastic tunneling of the Xe atom can appear only at certain resonant values 
of the bias voltage. At resonance the atom oscillates between a localized 
state on the surface (ipo), and a state close to tip (i^m), with a frequency 
attaining the maximum (~ l/2T max = 34.8 GHz) at the first resonance 
(U ~ —1.141 V). Thus, the exact treatment of the atom tunneling in the 
STM potential predicts an oscillatory behavior which cannot explain the 
irreversible switching mode. The observation of irreversible transfer with a 
constant rate indicates that in the present experiments the phase coherence is 
destroyed by the voltage pulse, or by other external factors, like dissipation. 
According to the quasi-classical results, the switching might be explained by 
a decoherence mechanism producing excitation and partial localization. 

The atomic QCO in STM might be observed as oscillations of the junction 
impedance, which changes when the atom moves from surface to the tip 
[11]. Therefore, accurate measurements of the resonant bias voltages and of 
the transfer time-scale could provide significant insight on the mechanism of 
decoherence and localization at the atomic level. 
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Figure Captions 

Fig. 1. The potential V at U = -1.141 V (solid), U = (dash), and the 
wave functions vpo, i^m (in units of A -1 / 2 ) as a function of distance. 
Fig. 2. The maximum value attained by p within 20 ps as a function of the 
bias voltage (A) and p as a function of time at the first resonance (B). 
Fig. 3. The polynomial potential V po i (solid) and the average potential V av 
(dash) for v = 0, 0.1, 0.5, 1, 1.5, 2, as a function of distance (A); xq C o (solid) 
and x c (dash) as a function of time (B). 

Fig. 4. The switching mode: phase-space orbit of the CM motion (A) and 
the CM coordinate as a function of time (B). 
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